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1 Introduction 

Some two hundred years ago, Joseph-Louis Lagrange (1736-1813) and Simeon De- 
nis Poisson (1781-1840) published articles which contained the first appearance of 
symplectic and Poisson structures, and of related conceptsQ The word symplectic, 
used for the first time with its modern mathematical meaning by Hermann Weyl (1885- 
1955) in his book The classical groups, first published in 1939 |29|, derives from a 
Greek word meaning complex. Weyl used it because the word complex, whose origin 
is Latin, had already a different meaning in mathematics. However, the concept of a 
symplectic structure is much older than the word symplectic since it appeared in the 
works of Lagrange, first in his 1808 paper I.15J about the slow variations of the orbital 
elements of the planets in the solar system, then again a few months later in [ 16 1 as a 
fundamental ingredient in the mathematical formulation of any problem in mechanics. 

Most modern textbooks present, as a first and fundamental example of a symplectic 
structure, the structure determined on the cotangent bundle of a smooth manifold by 
the exterior derivative of its canonical 1-form. It is in a slightly different context that 
the concept of a symplectic structure first appeared in the work of Lagrange, since it 



'l recall the memory of Nikolay Nekhoroshev, who passed away on 18 October 2008. He obtained many 
important results in the theory of dynamical systems, in particular some fundamental results on the stability in 
finite time intervals of almost completely integrable Hamiltonian systems, which may be seen as a continuation 
of the works of Lagrange and Poisson discussed in the present paper. 
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is on the manifold of motions of a mechanical system, rather than on the phase space 
of that system, i.e., the cotangent bundle of its configuration space, that he defined a 
symplectic structure. There are several reasons for considering Lagrange's point of 
view as more appropriate than the current one. For example, in most modem text- 
books, the conservation of the symplectic 2-form under the flow of a Hamiltonian 
vector field is presented as an important theorem in symplectic geometry, while this 
result had already been known to Lagrange, who considered it to be a direct conse- 
quence of the existence, on the manifold of motions of the system, of a well defined, 
time-independent symplectic structure. 

While Lagrange introduced the concept of a symplectic structure, Poisson defined 
the composition law today called the Poisson bracket. From June 1808 to February 
1810, five papers were published, two by Poisson and three by Lagrange, each paper 
improving on the results of the preceding one. 1 will present a reading of these works in 
the language of today's mathematicians, and I shall use modem notations and concepts 
when they can help to better understand Lagrange's and Poisson's ideas. My point of 
view will be that of a working mathematician, rather than that of a specialist in the 
history of mathematics. 

Section 2 describes two closely related concepts: the flow and the manifold of mo- 
tions of a smooth dynamical system. Section 3 describes the problem of the quantita- 
tive determination of the motion of the planets in the solar system, which was the main 
motivation for the work on dynamics by Lagrange and Poisson. Section 4 presents 
Lagrange's and Poisson's works on the method of varying constants. Finally, Section 
5 offers a modern account of this method, using today's notations and concepts. 



2 Flow and manifold of motions 

Let us consider a smooth dynamical system, i.e., an ordinary differential equation, on 
a smooth manifold M, 

On the right-hand side, X : i2 — > TM is a smooth, i.e., C°°, vector field, defined on an 
open subset H of M x M. The map (p, defined on an open interval of R, with values in 
M, is said to be a solution of this equation. For greater generality, we consider that, for 
t gR and X e M such that {t,x) €Q.,X{t,x), which is an element in TxM, may depend 
on t. We then say that X is a time-dependent vector field. Of course, a smooth vector 
field defined on M in the usual sense can be considered as a time-dependent vector 
field defined on K x M which, for each x G M, is constant on the line M. x {x}. 

2.1 The flow of a smooth differential equation 

The flow of a differential equation is the map, 

{t,to,Xo)>-^^{t,to,Xo), 

defined on a subset D of M x R x M, with values in M, such that, for each pair {to , xq ) G 
Q., the map t ^{t , to, xq) is the maximal solution of the differential equation which 
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takes the value xq at ? = ?o- Thus, by definition, 

d0{t,jo,xo) =x{t,^{t,to,xo)) , ^{to,to,xo) =xo. 

We recall that a solution of a differential equation, defined on an open interval /, is 
said to be maximal if it is not the restriction of a solution defined on an open interval 
strictly larger than /. 

It is well known that the map <I> is smooth; the subset D on which it is defined is 
open in M X R X M; for each {to,xo) € H, the set of reals t such that {t,to,xo) belongs 
to D is an open interval I{to,xo) which contains to- In addition, for any {to,xo) G H, fi 
and ?2 e M, 

^{t2,hMti,to,Xo)) = '^it2,to,Xo) ■ 

More precisely, if the left-hand side of this equaUty is defined, i.e., if (fi,fo,xo) £ D 
and (t2,t\,'P{ti,to,xo)) G D, then the right-hand side is defined, i.e., (f2,fo,xo) £ D, 
and the equality holds. Conversely, if the right-hand side and <I>(fi,fo,xo) are defined, 
i.e., if both (f2,fo,xo) &D and {ti,to,xo) G D, then the left-hand side is defined, i.e., 
(t2,ti,'i>{ti,to,xo)) G D, and the equality holds. 

The proof of the existence of maximal solutions, and therefore the proof of the 
existence of the flow of a smooth differential equation, rests on the axiom of choice. 



2.2 The manifold of motions of a smooth dynamical system 

The concept of manifold of motions of a smooth dynamical system is closely related 
to the concept of flow. The manifold of motions of a smooth dynamical system is the 
set M of all themaximal solutions t^<p{t) of the corresponding differential equation. 

As a set, M is the quotient of the open subset £2 of R x M, on which the time- 
dependent vector field X is defined, by the equivalence relation, 

(?2,X2) and (?i,xi) G SI are equivalent if {t2-t\.xi) belongs to the open subset!) of 
R X R X M on which the flow 4> is defined and 

X2 =4>(f2,fl,Xi). 

A smooth manifold structure on M can be defined as follows. An element uq&M 
is an equivalence class for the above-defined equivalence relation. Let (?o,xo) G i2 be 
an element of this equivalence class. According to the theorem of global existence and 
uniqueness of the maximal solution of a smooth ordinary differential equation which 
satisfies given Cauchy data, there exists an open neighborhood U{to.xo) of xq in M such 
that, for each x G there exists a unique maximal solution a which takes the 

value X att — to. We can use the map x i— > a to build a chart of M, whose domain is 
diffeomorphic to t/((Q,xo) • Using this construction for all {to,xo) G Q., we obtain an atlas 
of M, therefore a smooth manifold structure on this set, for which each point in M has 
an open neighborhood diffeomorphic to an open subset of M. The resulting smooth 
manifold structure of M need not be Hausdorff. 

When the flow <t> is defined on R x R x M, the manifold of motions M is globally 
diffeomorphic to M. But there is no canonical diffeomorphism of M onto M: the choice 
of a particular time ?o G R determines a diffeomorphism of M onto M which associates 
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with each motion a G M the point a{to) £ M. Of course this diffeomorphism depends 
on tQ. 

2.3 The modified flow 

The value <P{t,tQ,xo) of the flow <I> at time f depends on f and on the equivalence 
class fl e M of (fo,xo) G IR X M, rather than on fo and xq separately. Indeed, if (fc^co) 
and (fijjci) are equivalent, xi = <I>(fi ,fo,xo) and <I>(f,fi,xi) ~ <!>(?, fi,<I>(fi,fo,ji[:o)) = 
4>(f,fo,xo). The modified flow of the differential equation is the map, defined on an 
open subset of M x M, 

(f,fl) H^4>(r,fl) =4>(f,fo,Ji:o), 

where {to^xo) G M x M is any element of the equivalence class a <E M. As explained 
below, Lagrange used this concept in a paper [16,1 published in 1809. 

3 The motion of the planets in the solar system 

During the eighteenth century, the quantitative description of the motion of the planets 
in the solar system was a major challenge for mathematicians and astronomers. Let us 
briefly indicate the state of the art on this subject before the publications of Lagrange 
and Poisson. 

3.1 Kepler's laws 

In 1609, using the measurements obtained by the Danish astronomer Tycho Brahe 
(1546-1601), the German astronomer and physicist Johannes Kepler (1571-1630) dis- 
covered two of the three laws which very accurately describe the motion of the planets 
in the solar system. His discovery of the third law of motion followed in 1611. 

Kepler's first law 

As a first (and very good) approximation, the orbit of each planet in the solar system 
is an ellipse, with the sun at one of its foci. 

Kepler's second law 

As a function of time, the motion of each planet is such that the area swept by the line 
segment which joins the planet to the sun increases linearly with time. 

Kepler's third law 

The ratio of the squares of the rotation periods of two planets in the solar system is 
equal to the ratio of the cubes of the major axes of their orbits. 
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3.2 Orbital elements of a planet 



In Kepler's approximation, the motion of each planet in the solar system is a solution of 
a smooth dynamical system known as Kepler's problem, the motion of a massive point 
in a central attractive gravitational field, in Euclidean three-dimensional space. This 
problem, first formulated in mathematical terms by Isaac Newton (1642-1727), was 
solved by him in 1679 and published in 1687 in his Philosophiae naturalis principia 
mathematica ||22J . Each possible motion is determined by six quantities, called the 
orbital elements of the planet. Altogether these six elements constitute a system of 
local coordinates on the manifold of motions of Kepler's problem, and therefore the 
dimension of this manifold is 6. Let us explain why. 

First, we must determine the plane containing the planet's orbit. It is a plane con- 
taining the attractive center the sun. Such a plane may be determined, for example, 
by a unit vector with the sun as its origin, normal to this plane, i.e., a point of the 2- 
dimensional unit sphere centered at the sun. For this, we need 2 orbital elements. The 
choice of such a vector simultaneously determines an orientation of the orbital plane, 
which will be assumed such that the planet rotates counter-clockwise around the sun. 

We need two more orbital elements, which determine the orbit's shape and position 
in its plane. We know that this orbit is an ellipse with the sun as a focus. So the 
shape and position of the orbit are completely determined by the excentricity vector, 
discovered by Jakob Hermann (1678-1753), sometimes improperly called the Laplace 
vector or the Lenz vector (see A. Guichardet's paper |4|). It is the dimensionalless 
vector contained in the orbit's plane, directed from the attractive center towards the 
planet's perihelion, whose length is equal to the orbit's excentricity. 

We still need an orbital element to determine the size of this elliptic orbit. For 
example, we may choose the length of its major axis. 

Up to now we have seen that five orbital elements are needed to determine the 
planet's orbit. A sixth and final orbital element will determine the planet's position on 
its orbit. We may, for example, choose the point on this orbit at which the planet is 
at a fixed particular time. The second and third Kepler laws then fully determine the 
planet's position at all times, past, present and future. 

3.3 The manifold of motions of a planet in Kepler's approximation 

A modern and general description of the manifold of motions of a planet in Kepler's 
approximation was made by J.-M. Souriau [27|. He considered all possible motions, 
parabolic and hyperbolic as well as elliptic. By using a transform called regularization 
of collisions, he even included singular motions, in which the planet moves along a 
straight line until it collides, at a finite time, with the sun, and their analogues when 
time is inverted, in which a planet is, at a given time, ejected by the sun. This manifold 
is 6-dimensional. Due to the singular motions, it is non-Hausdorff. Other modern 
treatments of Kepler's dynamical system may be found in the books by V. Guillemin 
and S. Sternberg fS) and by B. Cordani |3|. 

Lagrange and Poisson were only interested in the elliptic motions of planets, not 
in parabolic or hyperbolic motions, which would be motions of comets rather than of 
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planets. In Kepler's approximation, the set of all elliptic motions of a planet, excluding 
singular motions, is an open, connected, Hausdorff submanifold of the manifold of all 
motions. We have seen above that the six orbital elements of the planet constitute a 
system of local coordinates on this manifold. If we choose a particular value of the 
time and consider the three coordinates of the planet and the three components of its 
velocity at that time, in any space reference frame, we obtain another system of local 
coordinates on the manifold of motions, and another way of showing that its dimension 
is 6. 

3.4 Beyond Kepler's approximation 

Kepler's approximation is only valid under the assumptions that each planet interacts 
gravitationally exclusively with the sun, and that its mass is negligible compared with 
that of the sun. In fact, even if one does not take into account the gravitational inter- 
action between planets, unless one assumes that the mass of the planet is negligible 
compared with that of the sun, its orbit is an ellipse whose focus is the center of mass 
of the system planet-sun, not the center of the sun. This center of mass is different for 
each planet. Therefore the planets have two kinds of gravitational interactions: their 
direct mutual interactions, and the interaction that each of them exerts on all the other 
planets through its interaction with the sun. 

To go beyond Kepler's approximation, astronomers and mathematicians used a 
very natural idea: each planet was considered to be moving around the sun on an 
ellipse whose orbital elements slowly vary in time, instead of remaining rigorously 
constant, as in Kepler's approximation. While astronomers increased the accuracy of 
their observations, from which they deduced tables for these slow variations, mathe- 
maticians sought to calculate them, using Newton's law of gravitational interaction. 
Let us briefly indicate some important stages of this search, which finally led to the 
mathematical discoveries made by Lagrange and Poisson. 

In 1773, Pierre-Simon Laplace (1749-1827) proved that, up to the first order, the 
gravitational interactions between the planets cannot produce secular variations in their 
periods, nor in the length of their orbit's major axis 1 191. Then in 1774, after reading 
Lagrange's paper discussed below, he calculated the slow variations of other orbital 
elements, the excentricity and the aphelion's position 1201 . In ||2TI he improved on his 
results for the planets Jupiter and Saturn. 

In 1774, Lagrange calculated the variations of the position of the nodes and or- 
bital inclinations of the planets UJJ. Then, in several papers presented to the Berlin 
Academy of Sciences between 1776 and 1784, he improved on Laplace's results, and 
determined the slow variations of other orbital elements. He distinguished between 
secular variations, non-periodic, or periodic with very long periods, which may be- 
come large with time l[l2][l3l, and periodic variations, which remain bounded for all 
times {TA\. He proved, with fewer approximations than Laplace, that the gravitational 
interactions between the planets cannot produce secular variations of their periods. 

Then, it seems that for more than 20 years, Lagrange ceased being interested in the 
subject, and he published no important new results on the motion of the planets. 
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On 20 June 1808, Poisson presented a paper, Sur les inegalites seculaires des 
moyens mouvements des planetes, to the French Academy of Sciences |25|, in which 
he removed a simpHfying assumption that had been made by Lagrange in his papers of 
the years 1776-1784 on the variation of the periods of the planets. 

Stimulated by Poisson's contribution, Lagrange returned to the problem in his 
Memoire sur la theorie des variations des elements des planetes, presented to the 
French Academy of Sciences on 22 August 1808 1 15|. We quote a passage from his 
introduction which shows that he clearly understood that Poisson's result was due to a 
still hidden mathematical structure. After recalling Laplace's important result of 1773, 
and the improvements he had obtained in 1776, he wrote. 

On n 'avait pas ete plus loin sur ce point; mais M. Poisson y a fait un pas de plus 
dans le Memoire qu 'il a lu il y a deux mois a la Classe, sur les inegalites seculaires des 
moyens mouvements des planetes, et dont nous avons fait le rapport dans la derniere 
seance. II a pousse 1' approximation de la meme formule jusqu'aux termes affectes des 
Carres et des produits des masses, en ayant egard dans cette formule a la variation des 

elements que j' avals regardes comme constants dans la premiere approximation 

il parvient d'une maniere ingenieuse a faire voir que ces sortes de termes ne peu- 
vent non plus produire dans le grand axe de variations proportionnelles au temps. 
...II me parut que le resultat qu 'il venait de trouver par le moyen des fonnules 
qui representent le mouvement elliptique etait un resultat analytique dependant de 
la forme des equations differentielles et des conditions de la variabilite des constantes, 
et qu'on devaity arrive r par la seule force de 1' Analyse, sans connaitre les expressions 
particulieres des quantites relatives a I'orbite elliptique. 

On 13 March 1809, Lagrange extended his method in the Memoire sur la theorie 
generale de la variation des constantes arbitraires dans tous les problemes de mecanique 

kMJ. 

On 16 October 1809, Poisson presented his paper, Sur la variation des constantes 
arbitraires dans les questions de mecanique ll26l . It is in this work, which is devoted 
to the subject that had been considered by Lagrange just a few months earlier, that he 
solved a question which had been left in abeyance by Lagrange, and introduced the 
composition law today called the Poisson bracket. 

On 19 February 1810, Lagrange presented his Second memoire sur la theorie de 
la variation des constantes arbitraires dans les problemes de mecanique fT7|. Here, 
he recognizes Poisson's contribution, but claims that the main ideas were already con- 
tained in his own previous paper. Lagrange included a simplified presentation of the 
method of varying constants in the second edition, published in 1 8 1 1 , of his Mecanique 
analytique ifTSll . 

It appears that, more than twenty years later, these works were still not fully un- 
derstood, since Augustin Louis Cauchy (1789-1857) had to give a very clear presen- 
tation of Lagrange's method in his Note sur la variation des constantes arbitraires 
dans les problemes de mecanique |2l. Published in 1837, this paper is a summary of 
a longer work he had presented in 1831 to the Academy of Sciences of Turin. In this 
work, Cauchy makes essential use of the formalism recently introduced by Sk William 
Rowan Hamilton (1805-1865) in his papers On a general method in Dynamics ^ and 
Second essay on a general method in Dynamics IT). 
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3.5 Remarks on the stability of the solar system 

At the beginning of the nineteenth century, the results obtained by Laplace, Lagrange 
and Poisson concerning the absence of secular variations of the major axis of orbits 
of the planets were regarded as a proof of the stability of the solar system. Today we 
know that these results were not rigorously established 1 1 1. In his famous article 123)1 . 
Jules Henri Poincare (1854-1912) proved that the problem is much more subtle than 
Laplace, Lagrange and Poisson had imagined. In |24|, vol. Ill, chapter XXVI, page 
140, he remarks that Poisson's results do not exclude the existence of terms of the form 
Af sin(ar + j3) in the expression for the variations of the major axis of a planet's orbit, 
f being the time. A, a and j3 denoting constants. Such terms can take very large values 
for large values of the time, although they vanish periodically. In fact, the problem of 
the stability of the solar system gave rise to modern research by many mathematicians, 
notably Andrei Nicolaievich Kolmogorov (1909-1987), Jurgen Moser (1928-1999), 
Vladimir Igorevich ArnoFd (1937), and Nikolai Nikolaievich Nekhoroshev (1946- 
2008). 

4 The method of varying constants 
4.1 Lagrange's paper of 1809 

Lagrange considers a mechanical system with kinetic energy, 

T = T{r,s, M, . . . ,r' ,s' ,u' . . .) , 

where r, s, u, . . . are independent real variables which describe the system's position in 
space. For a planet moving around the sun, these variables are the three coordinates of 
the planet in some reference frame. Let n be the number of these variables. In modern 
terms, n is the dimension of the configuration manifold. The quantities r', s', m', . . ., 
are the derivatives of r, s,u, . . ., with respect to the time, f , 

, dr I ds I du 

f — ^ S - — - ^ W — — ^ ^ • • • 

dt dt dt 

As a first approximation, Lagrange assumes that the forces which act on the system 
derive from a potential V, which depends on r, s, u, . . ., but not on the time derivatives, 
r', s', m', ... For a planet's motion, V is the gravitational potential due to the sun's 
attraction. The equations which determine the motion, established by Lagrange in his 
Mecanique analytique |,18J , are 

di [j? J ^ ~d7 ^ ~d7 ^ ' 

and similar equations in which r and r' are replaced by s and s', u and u', ... 

The general solution of this system of n second-order equations depends on the 
time f and on 2n integration constants. Lagrange denotes these constants by a, b, c, f, 
g, h,. . ., and writes this general solution as 

r^ r(t,a,b,c,f,g,h,...), s = s{t,a,b,c,f,g,h, . . .), u = ---. 
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For a planet's motion, the 2n integration constants a, b, c, f,g,h, ... are the orbital 
elements of the planet. 

At a second approximation, Lagrange assumes that the potential V does not fully 
describe the forces which act on the system, and should be replaced by V — il, where 
n may depend on r, s, u, and on the time t. For a planet's motion, Q. describes 
the gravitational interactions between the planet under consideration and all the other 
planets, which had been considered to be negligible in the first approximation. Q. 
depends on the time, because the planets which are the source of these gravitational 
interactions are in motion. The equations become 

d fdT\ dT dV _ da 

and similar equations in which r and r' are replaced by s and s', u and u', ... 
Lagrange writes the solution of this new system under the form 

r = r(t,a{t),b{t),c{t),f{t),g{t),h{t),...^ , 

and similar expressions for s, u, .... The function 

{t,a,b,c,f,g,h,...)^ r{t,a,b,c,f,g,h...) 

which appears in this expression, and the analogous functions which appear in the 
expressions of s, m, . . . are, of course, those previously found when solving the problem 
in its first approximation, with Q. set to 0. These functions are therefore considered as 
known. 

It only remains to find the 2n functions of the time f t-^ a{t), 1 1-^ b{t), These 

functions will depend on the time and on 2n arbitrary integration constants. 

4.2 Lagrange parentheses 

Lagrange obtains the differential equations which determine the time variations of 

these functions a{t), b{t), The calculations by which he obtains these equations 

are at first very complicated, and he makes two successive improvements, first in an 
Addition, then in a Supplement to his initial paper He finds a remarkable property: 
these equations become very simple when they are expressed in terms of quantities 
that he denotes by {a,b), (a,c), {a,f), {b,c), {b,f), .... Today, these quantities are 
still in use and they are called Lagrange parentheses. 

Lagrange parentheses are functions of a, b, c, f, g, h, They do not depend 

on time, nor on the additional forces which act on the system when Q. is taken into 
account. Jean-Marie Souriau ll27l |28| has shown that they are the components of the 
canonical symplectic 2-form on the manifold of motions of the mechanical system, in 

the chart of this manifold whose local coordinates are a, b, c, f, g, h, So Lagrange 

discovered the notion of a symplectic structure more than 100 years before that notion 
was so named by Hermann Weyl Ii29il . 
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We stress the fact that the Lagrange parentheses are relative to the mechanical 
system with kinetic energy T and applied forces described by the potential V . The 
additional forces described by Q. play no part in Lagrange's parentheses: the consid- 
eration of these additional forces permitted the discovery of a structure in which they 
play no part! 

At first, Lagrange obtained very complicated expressions for the parentheses (a, b), 
(a,c), {b^c), .... In the Aiic//f/o« to his paper (Section 26 of [ 17|), he obtained the much 
simpler expressions: 

( b) ^P'' drdpr ^dsdps dsdps^dudpu dudpu^ 

' da db db da da db db da da db db da ' 

and similar expressions for (fl,c), {b,c), . . . We have used the notations introduced by 
Hamilton \^J} and Cauchy |2| thirty years later, 

_dT _dT _dT 

Pr — 3-7 ; Ps — 3-7 , Pu — 3—7 , 

dr ds du 

while Lagrange used the less convenient notations, T' , T" and T'" , instead of pr, Ps 
and Pu. 

We recall that r, s, 11, . . . are local coordinates on the configuration manifold of the 
system, and r', s' , u' their partial derivatives with respect to time. The kinetic energy T, 
which depends on r, s, u, .. ., r', s', u' , . . ., is a function defined on the tangent bundle 
of the configuration manifold, which is called the manifold of kinematic states of the 
system. 

The map 

(r,i,M, . . .,r',s',u', . . .) 1-^ {r,s,u, . . . ,pr,Ps,Pu, ■ ■ ■) , 

called the Legendre transformation, is defined on the tangent bundle of the configu- 
ration manifold, and takes values in the cotangent bundle of this manifold, called the 
phase space of the system. When the kinetic energy is a positive definite quadratic 
form, this map is a diffeomorphism. This occurs very often, for example in the motion 
of a planet around the sun, the mechanical system considered by Lagrange. 

Since the integration constants a, b, c, f,g,h,... constitute a system of local 
coordinates on the manifold of motions, they completely determine the motion of the 
system. We again stress that this system is a first approximation, where Q. is set to 0. 
Therefore, for each time t, the instantaneous values of the quantities r, s, u, .. ., r', s', 
u', . . ., are determined as soon as a, b, c, f, g, h, ... are given. 

Conversely, the existence and uniqueness theorem for solutions of ordinary differ- 
ential equations (implicitly considered as obvious by Lagrange, at least for Kepler's 
problem whose solutions are explicitly known) shows that when the values of r, s, u, 
. . ., r' , s' , u', . . . at any given time f are known, then the motion is determined, so a, b, 
c, f,g,h, ... are known. 

In short, for each time f, the map which associates to a motion of coordinates 
(a, b, c, f, g, h,...) the values at time t of (r, s, u,...,r', s', u',...) is a diffeo- 
morphism from the manifold of motions onto the manifold of kinematic states of the 
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system. The composition of this diffeomorphism with the Legendre transformation 
yields, for each time t, a diffeomorphism from the manifold of motions onto the phase 
space, 

(a, b, c, f, g,h,...)\-^ {r{t), s{t), u{t), pr{t), Ps{t), Pu{t), ...), 

where r(f), s(f), u{t), Pr{t), Ps{t), Pu{t) are the values taken at time t by the corre- 
sponding quantities. 

The partial derivatives in the expression of the Lagrange parentheses are the partial 
derivatives of the diffeomorphism 

{a, b, c, f, g, h,...)h^ {r{t), s{t), u{t),. . . , pr{t), Ps{t), Pu{t), ■ ■ ■) 

where t is any value of the time, considered as fixed. 

Important remark 

The Lagrange parentheses {a,b) are defined when a complete system of local coordi- 
nates (a, b, c, f, g, h,.. .) has been chosen on the manifold of motions: the value of 
{a,b) depends not only on the values of the functions a and b on that manifold, but 
also on all the other coordinate functions, c, f, g, h, .... 

4.3 The canonical symplectic form 

Let us again consider the diffeomorphism 

(a, b, c, f, g,h,...)\-^ {r{t), s{t), u{t), pr{t), Ps{t), Pu{t), ...), 

where t is any value of the time, considered as fixed. The exterior calculus of differ- 
ential forms, created by Elie Cartan at the beginning of the twentieth century, did not 
exist in Lagrange's times. Today, with this very efficient tool, it is easy to prove that the 
Lagrange parentheses are the components of the pull-back by this diffeomorphism, on 
the manifold of motions, of the canonical symplectic 2-form of the cotangent bundle 
of the configuration manifold. In fact, 

{a,b)daAdb+ {a,c)daAdc-\ \- {b,c)db Adc -\ 

= f^da + ^db + • • • I A ( ^da + ^p-dh + • • • | 
\da db J \ da db J 

= dr Adpr + ds Adps + du Adpu-\ • 

The last expression is the well known formula for the components of a symplectic 
2-form in Darboux coordinates. 
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Another important remark 

Lagrange proved that, although they are defined by means of a diffeomorphism which 
depends on time, the parentheses he introduced do not depend explicitly on time: they 
are functions on the manifold of motions. When proving this result, Lagrange proved 
that the canonical symplectic 2-form on phase space is invariant under the flow of the 
evolution vector field on this space. 

4.4 Formulae for the variation of constants 

Lagrange proved that the derivatives with respect to time of the "constants that are 
varied" , a, b, satisfy 

,daj do. 

where, for short, 1 have written a,, 1 < / < 2n, instead of a, b, c, .. ., and where 1 have 
taken into account the skew-symmetry, {aj,ai) — —{ai.aj). 

Lagrange indicates that by solving this linear system, one obtains something like 

dat 2« 
dt oaj 

He explains that the Lij are functions of the a,- which do not depend explicitly on 
time. In modern terms, the Ljj are functions defined on the manifold of motions. But 
Lagrange does not state their explicit expressions. That would be done by Poisson a 
few months later 



4.5 Poisson's paper of 1809 and the Poisson bracket 

When he was a student at the Ecole Polytechnique, Poisson attended lectures by La- 
grange. In a paper read before the French Academy of Sciences on 16 October 1809 
||26| . he added an important ingredient to Lagrange's method of varying constants. He 
introduced new quantities, defined on the manifold of motions, which he denoted by 
{a,b), {a,c), . . . These quantities are not the Lagrange parentheses. Today, they are 
called the Poisson brackets. In his paper, Poisson also uses Lagrange parentheses but 
he denotes them differently, by [a,b] instead of {a,b), [a,c] instead of (a,c), etc. 

We shall retain Lagrange's notations {a,b), {a,c), . . . for the Lagrange parentheses 
and we will denote the Poisson brackets by {a,b}, {a,c}, 

The expression of the Poisson brackets is 

da db da db ^ da db da db ^ da db da db ^ 
' dprdr dr dp,- dps ds dsdps dpudu dudpu 

Of course {a,c}, {fe,c}, ... are given by similar formulae. We observe that in these 
formulae there appear the partial derivatives of the local coordinates a, b, c, .. . on the 
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manifold of motions, considered as functions of the dynamical state of the system at 
a fixed time, t. The independent variables which describe this dynamical state are the 
values, at time t, of the quantities r, pr, s, ps, u, /?„, 

The above formula is the well known expression of the Poisson bracket of two 
functions a and b defined on a symplectic manifold, in Darboux coordinates. 



4.6 Poisson brackets versus Lagrange parentheses 

Let us compare the Poisson bracket 

^ da db da db ^ da db da db ^ da db da db ^ 

' dpr dr drdpr dps ds ds dps dpi,du dudpi, 

with the Lagrange parenthesis 

^ dr dp,- drdpr^dsdps dsdp^^dudpu dudpu^ 
' da db db da da db db da da db db da 

We see that these formulae involve the partial derivatives of two diffeomorphisms 
which are inverses of one another: the Poisson bracket involves the partial derivatives 
of the coordinates a,b, ... on the manifold of motions with respect to the coordinates 
r, i, r', i', ... on the phase space, while the Lagrange parenthesis involves the partial 
derivatives of r, s, /, s' , . . . with respect to a,b, . . . 

We have seen that Lagrange's parentheses {a,b), {a,c), are the components 
of the symplectic 2-form on the manifold of motions, in the chart of this manifold 

whose local coordinates are a, b, c, The Poisson brackets {a,b}, {a,c}, . . ., are the 

components in the same chart of the associated Poisson bivector 

The matrix whose components are the Lagrange parentheses {a,b), {a,c), . . ., and 
the matrix whose components are the Poisson brackets {a,b}, {a,c}, . . ., are inverses 
of one another. This property was clearly stated by Cauchy in his paper [2], read 
before the Academy of Turin on 11 October 1831, 22 years after the publication of 
Lagrange's and Poisson's papers. 



4.7 Remark about the Poisson theorem 

The Poisson theorem states that the Poisson bracket of two first integrals, i.e., functions 
which remain constant on each trajectory of the system, is a first integral. Today, 
this result is often presented as a consequence of the Jacobi identity. This identity 
was not known to Lagrange, nor to Poisson, who considered the constancy of the 
Poisson bracket of two first integrals to be due to the fact that it is a function defined 
on the manifold of motions. The Poisson bracket can indeed be defined for any pair of 
smooth functions on the manifold of motions, and it is still a function defined on that 
manifold which depends only on the two given functions. As stated above, contrary 
to the Poisson bracket, the Lagrange parenthesis can be defined only for a pair of 
functions which are part of a complete system of coordinate functions, and not for a 
pair of smooth functions in general. 
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4.8 The Jacobi identity 



Lagrange and Poisson observed the skew-symmetry of their parentheses and brackets, 
but said nothing about the Jacobi identity for the Poisson brackets, nor about the rela- 
tions between Lagrange's parentheses expressing the fact that they are the components 
of a closed 2-form. 

Discovered by the German mathematician Carl Gustav Jacob Jacobi (1804-1851) 
(HI [To I, the identity that bears his name involves three arbirary smooth functions /, g 
and h defined on a symplectic (or a Poisson) manifold, 

{/, {gM} + {g. {hj}} + {K {f.g}} =0, 

or three smooth vector fields X, Y and Z defined on a smooth manifold, 

[X,[Y,Z]] + [Y,[Z,X]] + [Z,[X,Y]\ =0. 

Jacobi understood the importance of this identity, which later played an important part 
in the theory of Lie groups and Lie algebras developed by the Norwegian mathemati- 
cian Marius Sophus Lie (1842-1899). 



4.9 Lagrange's paper of 1810 

In his paper ifTTll . Lagrange gives simpler expressions of his previous results, using 
Poisson brackets. He writes the differential equations which determine the time varia- 
tions of the "constants" a, b, . . ., in the form 

—— = > \ai,aj \— — . \ < i<2n. 
dt ^, ^ ^' da; - - 

Here 1 have denoted the constants by a,, 1 < / < 2n, a notation that allows the use of 

2n 

the symbol ^ for a more concise expression. Lagrange used longer expressions in 

i=l 

which the constants were denoted by a, b, c, f,g,h,... 

Let us observe that Lagrange could have written his equations in the simpler form, 

— - = , 1 < / < 2n , 

dt 

since Q. can be considered as a function defined on the product of the manifold of 
motions with the factor K, for the time. Therefore the Poisson bracket {fl,,i2} can be 
unambiguously defined as 

{fl,-,n} = y{ai,aj}—. 

. , da ; 

Lagrange did not use this simpler expression, nor did Poisson in his paper of 1809. 
Both used the Poisson bracket only for coordinate functions a,, not for more general 
functions such as CI. 
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4.10 Cauchy's paper of 1837 



This short paper of 6 pages, pubHshed in the Journal de Mathematiques pares et ap- 
pUquees, is extracted from the longer paper presented by Cauchy before the Academy 
of Turin on 11 October 1831. Its title is almost the same as those of the papers by 
Lagrange and Poisson. 

Cauchy very clearly explains the main results due to Lagrange and Poisson, using 
Hamilton's formalism. However, he does not write Poisson brackets with the function 
n (which is denoted by R in his paper). 

Cauchy proves, without using the word matrix, that the matrix whose coefficients 
are Lagrange's parentheses of some coordinate functions and the matrix whose coef- 
ficients are the Poisson brackets of the same coordinate functions are inverses of one 
another. 



5 Varying constants revisited 

I shall now present, in modern language and with today's notations, the main results 
due to Lagrange and Poisson concerning the method of varying constants. 1 will follow 
Cauchy's paper of 1837. 

Let (M, 03) be a 2n-dimensional symplectic manifold, with a Hamiltonian function 
Q : M X M — > R which may be time-dependent {Q is the notation used by Cauchy). Let 
Mo be the manifold of motions of this Hamiltonian system, and let <I> : M x Mq M, 
{t,a) I— > a) be the modified flow of the Hamitonian vector field associated with Q 
(in the sense of Section l23T l. The easiest way of writing Hamilton's equation is the 
following. For each smooth function, g : M ^ M, 

" ^{Q,g}{^{t,a)). 

We now assume that the system's true Hamiltonian is Q+R instead of Q, where R 
may be time-dependent. 

The aim of the method of varying constants is to transform the modified flow of the 
Hamiltonian vector field associated with Q into the modified flow of the Hamiltonian 
vector field associated with Q + R. More precisely, the aim is to find a map ^ : E x 
Ml Mq, {t,b) t-^ a = ^'{t,b), where Mi is the manifold of motions of the system 
with Hamiltonian Q + R, such that {t,b) ^ <P(t ,^'{t ,b)) is the modified flow of the 
vector field with Hamiltonian Q + R. 



These maps must satisfy, for any smooth function g : M 

(go$(f,>P(f,/7))) = {Q + R,g}(^^t,^it,l 

For each value fo of the time f , 

d 



d 
dt 



d 
dt 



'='0 
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When to is fixed, {t,W{to,b)) ^ 4>(f,^(fo,^)) is the flow of the vector field with 
Hamiltonian Q. Thus the first term of the right-hand side is 

^(^o$(f,^(fo,/7))) [^^^={e,g}($(fo,^(fo,^) 

Therefore the second term of the right-hand side must be 

j^{go^to,nt,b))) [^^^ = [{Q+R,g}-{Q,g})(^{to,nto,b))) 

= {R,g}M{^to,'V{to,b))) 

= {7?0$,„,^o$,jMo(^(fO,^)). 

The Poisson bracket of functions on M is denoted by { , } when there is no risk of 
confusion, and by { , }m when we want to distinguish it from the Poisson bracket of 
functions defined on Mq, which is denoted by { , }mo- For the last equahty, we have 
used the fact that : A/q — * M is a Poisson map. 

The function go= 8° can be any smooth function on Mq, so the last equality 
may be written as 



= {Ro<P,^,go}Mo{'i'{to,b)). 

t=t(j 



dW{t,b)\ d{8o{'i'it,b) 

dgo, 



dt 



dt 



Now, to may take any value, so the last equation proves that for each b € Mi, the 
manifold of motions of the system with Hamiltonian Q + R, the map 1 1— > ^{t,b) is an 
integal curve, lying in the manifold Mq of motions of the system with Hamiltonian Q, 
of the Hamiltonian system with the time-dependent Hamiltonian, 

{t,a)l-^R[t,<^{t,a)) , {t,a)GRxMo. 

This is the result discovered by Lagrange circa 1808. It is an exact result, not an 
approximate one. However, when the method is used for the determination of the 
motion of a given planet of the solar system, the potential R depends upon the positions 
of all the other planets, which are not known exactly. Therefore the method must be 
used in conjunction with successive approximations, the value of R used at each step 
being that deduced from the calculations made at the preceding steps. 
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